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Abstract 

^»^ We consider the problem of minimizing L -^/^^ + K^{s) ds for a planar curve having fixed 

initial and final positions and directions. The total length £ is free. Here s is the variable of 

1^3 arclength paramctrization, K{s) is the curvature of the curve and ^ > a parameter. This 

y—i problem comes from a model of geometry of vision due to Pctitot, Citti and Sarti. 

QQ We study existence of local and global minimizcrs for this problem. We prove that if for a 

^Sj certain choice of boundary conditions there is no global minimizer, then there is neither a local 

minimizer nor a geodesic. 

We finally give properties of the set of boundary conditions for which there exists a solution 
to the problem. 
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1 Introduction 

In this paper we are interested in the following variational problem^: 



cn (P) Fix ^ > and {xin,yin,din)Axfin,yfin,Ofin) G M^ X S^ . On the space of (regular enough) 

^^ planar curves, parameterized by plane-arclength^ find the solutions of: 

OO 

O x(0) = (xj„,yj„), :si{£) = {xfin,yfin), 

cn 
en 

o 

(N 



x(0) = (cos(6'j„),sin(6'jn)), x(^) = (cos(6'/j„),sin(6'/j„)), 

ri 

/ \/S,'^ + K'^{s) ds — )• min (here £ is free.) (1) 

Jo 



.^ Here K = , .^^.^.^^^ is the geodesic curvature of the planar curve x(-) = (x (•),?/(•)). This problem 

rN comes from a model proposed by Petitot, Citti and Sarti (see [9, 21, 22, 26] and references therein) 

C^ for the mechanism of reconstruction of corrupted curves used by the visual cortex VI. The model 

is explained in detail in Section 2. 

It is convenient to formulate the problem (P) as a problem of optimal control, for which the 
functional spaces are also more naturally specified. 



^In this paper, by S*^ we mean R/ ~ where 6 -^ 6' ii 9 ^ 6' + 2mT, n G Z. By P^ wc mean E/ ^ where 9 ^ 9' li 
9 = 9' +nn,nGZ. 

^Here by plane-arclength we mean the arclength in K^, for which we use the variable s. Later on, we consider 
also parameterizations by arclength on R^ x S^ or M^ x P^ , that we call sub-Riemannian arclength (sR-arclength for 
short), for which we use the variable t. We will also use the variable t for a general parametrization. 



(Pcurve) Fix ^ > and {xin,yin,Oin),{xfin,Vfin,Gfin) G K^ X 5*^ III the space of integrable 
controls v{-) : [0,£] — )■ M, find the solutions of: 

(x(0),y(0),6l(0)) = {Xin,yin,Oin), {x{i),y{£),e{£)) = {x fin, y fin, fin), 

g(s) \ / corns)) \ / \ 

g(.) = sin(^(s)) ]+v{s){ (2) 

fsis)) \ J \lj 

/ V '^^ + K{s)'^ ds = \/S,'^ + f (s)^ ds — )■ min (here £ is free) (3) 

Jo ^0 

Since in this problem we are taking v{-) G L^{[0,i]), we have that the curve q{.) = {x{-),y{-),6{-)) : 
[0,£] -^ M? X S^ is absolutely continuous and the planar curve x(-) := {x{-),y{-)) : [0,i] — )■ M? is in 
1^2,1 ([o,£]). 

Remark 1 Notice that the function y .^^ + K"^ has the same asymptotic behaviour, for K ^ and 
for K ^ oo of the function (p{K) introduced by Mumford and Nitzberg in their functional for image 
segmentation (see [20]). 

The main issues we address in this paper are related to existence of minimizers for problem 
(Pcurve)- More precisely, for (Pcurve) the first question we are interested in is the following: 

Ql) Is it true that for every initial and final condition, the problem (Pcurve) admits a global 
minimum? 

In [6] it was shown that there are initial and final conditions for which (Pcurve) does not admit 
a minimizer. More precisely, it was shown that there exists a minimizing sequence for which the 
limit is a trajectory not satisfying the boundary conditions. See Figure 1. 
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Figure 1: Minimizing sequence g„ converging to a non-admissible curve q (angles at the begin- 
ning/end) . 



From the modeling point of view, the non-existence of global minimizers is not a crucial issue. 
It is very natural to assume that the visual cortex looks only for local minimizers, since it seems 
reasonable to expect that it primarly compares nearby trajectories. Hence, a second problem we 
address in this paper is the existence of local minimizers for the problem (Pcurve)- More precisely, 
we answer the following question: 

Q2) Is it true that for every initial and final condition the problem (Pcurve) admits a local 
minimum? If not, what is the set of boundary conditions for which a local minimizer exists? 

The main result of this paper is the following. 



Theorem 2 Fix an initial and a final condition qin = {xin, yin, din) o,nd qjin = (xfin, Hfin-, ^fin) *^ 
R^ X S*^ . The only two following cases are possible. 

1. There exists a solution (global minimizer) for (Pcurve) from qin to qfm- 

2. The problem (Pcurve) from q^ to qfm does not admit neither a global nor a local minimum 
nor a geodesic. 

Both cases occur, depending on the boundary conditions. 

We recall that a curve q{.) is a geodesic if for every sufficiently small interval [^1,^2] C Dom(q{.)), 
the curve Q'(-)|[ti,t2] i^ ^ minimizer between q{ti) and 9(^2)- 

One of the main interests of (Pcurve) is that it admits minimizers that are in W^'^ but are not 
Lipschitz, as we will show in Section 5.2. As a consequence, controls lie in L^ but not in L°°. This 
is an interesting phenomenon for control theory: indeed, to find minimizers, one usually applies 
the Pontryagin Maximum Principle (PMP in the following) , that is a generalization of the Euler- 
Lagrange condition. But the standard formulation of the PMP holds for L°° controls; this obliges 
us to use a generalization of the PMP for (Pcurve)) that we discuss in Section 5.1. Details of this 
interesting aspect of (Pcurve) are given in Section 5.2. This also explains the reason for which we 
need to define variational problems, global and local minimizers in the space VF^'^, see Section 3. 

The second sentence of Q2 is interesting, since one could compare the limit boundary conditions 
for which a mathematical reconstruction occurs with the limit boundary conditions for which a 
reconstruction in human perception experiments is observed. Indeed, it is well known from human 
perception experiments that the visual cortex VI does not connect all initial and final conditions, 
see e.g. [22]. With this goal, we have computed numerically the configurations for which a solution 
exists, see Figure 2. 

The structure of the paper is as follows. In Section 2 we briefiy describe the model by Petitot- 
Citti-Sarti for the visual cortex VI. We state it as a problem of optimal control (more precisely 
a sub-Riemannian problem), that we denote by (Pprojective)- The problem (Pcurve) is indeed 
a modified version of (Pprojective)- In Section 3 we recall definitions and main results in sub- 
Riemannian geometry, that is the main tool we use to prove our results. In Section 4 we define 
an auxiliary mechanical problem (crucial for our study), that we denote with (Pmec); and study 
the structure of geodesies for it. In Section 5 we describe in detail the relations between problems 
(Pcurve), (Pprojective) and (Pmec), with an emphasis on the connections between the minimizers 
of such problems. In Section 6 we prove the main results of the paper, i.e. Theorem 2. 

2 The model by Petitot-Citti-Sarti for VI 

In this section, we recall a model describing how the human visual cortex VI reconstructs curves 
which are partially hidden or corrupted. The goal is to explain the connection between reconstruc- 
tion of curves and the problem (Pcurve) studied in this paper. 

The model we present here was initially due to Petitot [21, 22], based on previous work by 
Hubel-Wiesel [17] and Hoffman [15], then refined by Citti and Sarti [9, 26], and by the authors of 
the present paper in [8, 11, 12]. It was also studied by Hladky and Pauls in [14]. 

In a simplified model^ (see [22, p. 79]), neurons of VI are grouped into orientation columns, 



■'For example, in this model we do not take into account the fact that the continuous space of stimuli is implemented 
via a discrete set of neurons. 
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Figure 2: Configurations for wliicfi we have existence of minimizers with ^ = 1. For other ^ 7^ 1, one 
can recover the corresponding figure via dilations, as explained in Remark 20. Due to invariance 
of the problem under rototranslations on the plane, one can always assume that qin = (0,0,0). 

4, with y > 0. The case y < can be 



Left: We study the cases x'l-^ 



+ yjin 



1 and x2 .„ 



+ yh 



recovered by symmetry. In the case x'i-^ + y'j^^ = 1 minimizing curves are also shown. Right: For 
each point on the right half-plane, we give the set of configurations for which we have existence of 
minimizers. 



each of them being sensitive to visual stimuli at a given point of the retina and for a given direction 
on it. The retina is modeled by the real plane, i.e. each point is represented by {x,y) G M?, 
while the directions at a given point are modeled by the projective line, i.e. 6 ^ P^. Hence, the 
primary visual cortex VI is modeled by the so called projective tangent bundle PTE? := E? x P^. 
From a neurological point of view, orientation columns are in turn grouped into hypercolumns, 
each of them being sensitive to stimuli at a given point (x, y) with any direction. In the same 
hypercolumn, relative to a point (x, y) of the plane, we also find neurons that are sensitive to 
other stimuli properties, like colors, displacement directions, etc... In this paper, we focus only 
on directions and therefore each hypercolumn is represented by a fiber P^ of the bundle PTM?. 
Orientation columns are connected between them in two different ways. The first kind is given by 
vertical connections, which connect orientation columns belonging to the same hypercolumn and 
sensible to similar directions. The second is given by the horizontal connections, which connect 
orientation columns in different (but not too far) hypercolumns and sensible to the same directions. 
See Figure 3. 
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Figure 3: A scheme of the primary visual cortex VI. 



In other words, when VI detects a (regular enough) planar curve (a;(-),y(-)) 
computes a "lifted curve" in PTM? by including a new variable 9{-) : [0,T] — ) 
PFi'i([0,r]), which satisfies: 



[0, T] -^ M2 it 
P^, defined in 






U[T 



cos(^(r)) 

sin(6'(T)) 





+ v{t) 




for some 



u,v 



[0,T] 



(4) 



The new variable 9{.) plays the role of the direction in P^ of the tangent vector to the curve. Here 
it is natural to require u{-),v{-) G L^([0,r]). This specifies also which regularity we need for the 
planar curve to be able to compute its lift: we need a curve in W'^'^. 

Observe that the lift is not unique in general: for example, in the case in which there exists 
an interval [Ti,r2] such that gf (t) = j^{t) = for all r G [Ti,r2], one has to choose n = on the 
interval, while the choice of v is not unique. Nevertheless, the lift is unique (modulo L^) in many 
relevant cases, e.g. if ^(t)^ + ^(t)^ = for a finite number of times r G [0,T]. 

In the follow^ing ^ve call a planar curve a liftable curve if it is in W"^'^ and its lift is 
unique. 

Consider now a liftable curve (x(-),y(-)) : [0,T] — t- M? which is interrupted in an interval 
(a, 6) C [0,r]. Let us call {xin,yin) ■= (x(a),y(a)) and {x fin, y fin) ■= {x{b),y{b)). Notice that 
the limits 9in := limT-i-a 6'(r) and 9 fin := limT-j^f, 6'(t) are well defined, since 9{.) is an absolutely 
continuous curve. In the model by Petitot, Citti, Sarti and the authors of the present article 



[8, 9, 23], the visual cortex reconstructs the curve by minimizing the energy necessary to activate 
orientation columns which are not activated by the curve itself. This is modeled by the minimization 
of the functional 

rb 

J = {C u{t) + v{t) ) (ir — 7- min, (here a and b are fixed). (5) 

J a 

Indeed, ^^u(r)^ (resp. f (t)^) represents the (infinitesimal) energy necessary to activate horizontal 
(resp. vertical) connections. The parameter ^ > is used to fix the relative weight of the horizontal 
and vertical connections, which have different phisical dimensions. The minimum is taken on the set 
of curves which are solution of (4) for some u{-),v{-) G L^Ha, b]) and satisfying boundary conditions 

{x{a),y{a),e{a)) = {xin,yin,Oin), {x{b),y{b),e{b)) = (x fin, y fin, Of in)- 

Minimization of the cost (5) is equivalent to the minimization of the cost (which is invariant by 
reparameterization) 

C= f V^Mr? + v{Ty dT= f \\±{T)We + K{Ty (It, 

J a J a 

where x = (x, y) and with b > a fixed. See a proof of such equivalence in [18]. 
We thus define the following problem: 

(Pprojective): Fix ^ > and ixin,yin,Oin), (x fin, y fin, fin) G M^ X P^ . In the space of integrable 
controls u{-),v{-) : [0,T] — )■ M, find the solutions of: 

(a;(0), y(0), 6^(0)) = (xj„, yj„, Bin), {x{T), y{T), e{T)) = (x/i„, y/i„, 6 fin), 

I cos(0(r)) \ 
= u{t) sin(0(r)) +v{t) 

C= [ ^C'^u{t)^ + t;(T)2 dr = I \\±{T)\\y/^^^TK{Tf dT ^ uAn (here T is free) 

Observe that here 6 ^ P^, i.e. angles are considered without orientation^. 

The optimal control problem (Pprojective) is well defined. Moreover, it is a sub-Riemmanian 
problem, see Section 3. We have remarked in [8] that a solution always exists. We have also studied 
a similar problem in [7], when we deal with curves on the sphere S'^ rather than on the plane M?. 

One the main interests of (Pprojective) is the possibility of associating to it a hypoelliptic 
diffusion equation which can be used to reconstruct images (and not just curves), and for contour 
enhancement. This point of view was developed in [8, 9, 11, 12]. 

However, its main drawback (at least for the problem of reconstruction of curves) is the existence 
of minimizers with cusps, see e.g. [6]. Roughly speaking, cusps are singular points in which velocity 
changes its sign. More formally, we say that a curve trajectory {q{.),(u{.),v{.))) has a cusp at 
r G [0, T] if u{t) changes its sign in a neighbourhood^ of f. Notice that in a neighborhood of a 
cusp point, the tangent direction (with no orientation) is well defined. A minimizer with cusps is 
represented in Figure 4. 





''Notice that the vector field (cos 6, sin 9, 0) is not continuous on PTM.'^. Indeed, a correct definition of (Pprojective) 
needs two charts, as explained in detail in [8, Remark 12]. In this paper, the use of two charts is implicit, since it 
plays no crucial role. 

^More precisely, it exists £ > such that u{a)u{b) < for almost every a G (r — e, r) , 6 G (r, f + e). 



Figure 4: A minimizer with two cusps. 



However, to our knowledge, the presence of cusps has not been observed in human perception 
experiments, see e.g. [22]. For this reason, people started looking for a way to require that no 
trajectories with cusps appear as solutions of the minimization problem. In [9, 12] the authors 
proposed to require trajectories parameterizated by spatial arclength, i.e. to impose ||x|| = u = 1. 
In this way cusps cannot appear. Notice that assuming u = 1, directions must be considered with 
orientation, since now the direction of x is defined in S^ . In fact, cusps are precisely the points 
where the spatial arclength parameterization breaks down. By fixing n = 1 we get the optimal 
control problem (Pcurve) on which this paper is focused. 

Remark 3 We also define an "angular cusp" as follows: we say that a pair trajectory-control 
{q{.), {u{-), f (.))) has an "angular cusp" atf£ [a, b] if there exist a neighbourhood B := (r — e, f +e) 
such that u{t) = on B and9{f — e) ^ 9{f + e). Angular cusps are of the kind q{T) = {xo,yo,9o + 
IoVicr)da). 

The minimum of the distance between (xo,yO)^o) o.'i^d {xo,yQ,6i) with arbitrary 9o,Oi is realized 
by such kind of trajectories. This is the only interesting case in which we need to deal with such 
trajectories. Indeed, even assuming that a solution q{.) of (Pprojective) satisfies ui = on a 
neighbourhood t only, then analyticity of the solution^ implies that ui = on the whole [0,T], and 
hence q{.) steers {xo,yo,6o) to some {xo,yo,9i). 

3 Optimal control 

In this section, we give the fundamental definitions and results from optimal control, and the 
particular cases of sub-Riemannian problems, that we will use in the following. For more details 
about sub-Riemannian geometry, see e.g. [5, 13, 19]. 



3.1 Minimizers, local minimizers, geodesies 

In this section, we give main definitions of optimal control. Observe that we deal with curves in 
the space W^'^ to deal with the problem (Pcurve); see Section 5.2. 



^Analyticity of the solution is proved below, see Remark 13. 



Definition 4 Let M be an n dimensional smooth manifold and fu'-Q*-^ fu{<l) G TqM he a 1- 
parameter family of smooth vector fields depending on the parameter u G M™". Let f^:Mx W^ — t- 
[0, +oo) be a smooth function of its arguments. We call variational problem (denoted by (VP) 
for short) the following optimal control problem 



T 



fu(T){q{T-)), (6) 



/*'(g'(r), u(t)) dr — )■ mill, T free (7) 


q{0)=qo, q{T)=qi, (8) 

u{-) G UT>o^i([0,r],M™), q{-) G Ut>o^^''H[0, T], M) (9) 

Following [28, Ch. 8], we endow VJT>QW^^^{[f),T],M) with a topology. 

Definition 5 Let qi{.),q2i-) e Ut>oW^^'^{[0,T], M), with qi defined on [0, Ti] and q2 on [0,T2]. 
Extend qi on the whole time-interval [0,max{Ti,T2}] by defining qi{t) := qi{Ti) for t > Ti, and 
similarly for q2- We define the distance between qi{) and (?2(-) os 

Iki(-) - 92(.)IIh^m := ki(0) - 92(0)1 + ||(?i(.) - g2(.)llLi- 

From now on, we endow Ut>oW^^'"'^([0, T], M) with the topology induced by this distance. It is clear 
that this distance is induced by the norm in H^^'^. For more details, see [28, Ch. 8]. 
We now give definitions of minimizers for (VP). 

Definition 6 We say that a pair trajectory-control {q{-),u{-)) is a minimizer if it is a solution of 

(VP). 

We say that it is a local minimizer if there exists an open neighborhood Bqt.\ ofq{) in Ur>oVF^'^([0, T], 

endowed with the topology defined above, such that all {q{-), u{-)) satisfying (6)-(8), with q{-) G Bqf\, 

have a larger or equal cost. 

We say that it is a geodesic if for every sufficiently small interval [^1,^2] C Dom{q{-)), the pair 



(g(-),ti(-))|[i^^t2] ^^ ^ minimizer of j^ /°(g(r), ii(T)) dr from q{ti) to 9(^2) with T fi 



ree. 



Remark 7 It is interesting to observe that, in general, one can have the same trajectory q{.) 
realized by two different controls ui{.),U2{-)- For this reason, one has to specifiy the control to have 
the cost of a trajectory. Nevertheless, for the problems studied in this article, it is easy to prove 
that, for a given trajectory q{.) G H^^'"'^([0,T],M™) satisfying (6) for some control u{.), then such 
control is unique. 

In this paper we are interested in studying problems that are particular cases of (VP), see 
Section 5. In particular, we study the problem (Pmec) defined in Section 4, that is a 3D contact 
problem (see the definition below). For such problem we apply a standard tool of optimal con- 
trol, namely the Pontryagin Maximum Principle (described in the next section), and then derive 
properties for (Pcurve) from the solution of (Pmec)- 

3.2 Sub-Riemannian manifolds 

In this section, we recall the definition of sub-Riemannian manifolds and some properties of the 
corresponding Carnot-Caratheodory distance. We recall that sub-Riemannian problems are special 
cases of optimal control problems. 



Definition 8 A sub-Riemannian manifold is a triple (M, A,g), where 

• M is a connected smooth manifold of dimension n; 

• A. is a Lie bracket generating smooth distribution of constant rank m < n; i.e., k is a smooth 
map that associates to q £ M an m-dim subspace A{q) of TqM , and \/ q £ M, we have 

span{[Xi, [. . . [Xk-i,Xk] . . .]]{q) \ Xi G Vec(M) and Xi{p) G k{p) V p G M} = T^M. 

Here Vec(M) denotes the set of smooth vector fields on M. 

• g,q is a Riemannian metric on A.{q), that is, smooth as a function of q. 

The Lie bracket generating condition (10) is also known as the Hormander 
condition, see [16]. 

Definition 9 A Lipschitz continuous curve q{.) : [0,T] -^ M is said to be horizontal if q{T) G 
A(g(r)) for almost every r G [0,T]. Given a horizontal curve q{.) : [0,T] — )• M, the length of q{.) is 



m-)) = j ^%,irMT)A{T))dT. (10) 



'0 
The distance induced by the sub-Riemannian structure on M is the function 

d{qo,qi) = iiif {/((7(.)) I 9(0) = qo, <l{T) = qi, q{.) horizontal Lipschitz continuous curve} . (11) 

Notice that the length of a curve is invariant by time-reparametrization of the curve itself. 

The hypothesis of connectedness of M and the Lie bracket generating assumption for the dis- 
tribution guarantee the finiteness and the continuity of d{-,-) with respect to the topology of M 
(Rashevsky-Chow's theorem; see, for instance, [4]). 

The function d{-, •) is called the Carnot-Caratheodory distance. It gives to M the structure of 
a metric space (see [5, 13]). 

Observe that (Pprojective) and (Pmec) defined in Section 4 are both sub-Riemannian problems. 
Indeed, defining 

/ cos6l \ / 

Xi= sin^ , X2= 

one has that (Pprojective) is sub-Riemannian with M = PTM?, kg = span{Xi(g), X2{q)} and g{q) 
such that Xi{q),X2{q) is an orthonormal basis. For (Pmec)i simply replace M = SE{2). 

3.3 The Pontryagin Maximum Principle on 3D contact manifolds 

In the following, we state some classical results from geometric control theory which hold for the 
3D contact case. For simplicity of notation, we only consider structures defined globally by a pair 
of vector fields, that are sometimes called "trivialized structures". 



Definition 10 (3D contact problem) Let M be a 3D manifold and let Xi,X2 be two smooth 
vector fields such that dim(Span{Xi, X2, [Xi, X2]}{q))=3 for every q G M. The variational problem 



T 

min (11) 



q = UiXi + U2X2, qiO) = qo, q{T) = qi, I yj {ui{t)Y + {u2{t)Y dT 

Jo 

is called a 3D contact problem. 

Observe that a 3D contact manifold is a particular case of a sub-Riemannian manifold, with ▲ = 
span{Xi,X2} and Sq(t) is uniquely determined by the condition ggr^\{Xi,Xj) = 6ij. In particular, 
each 3D contact manifold is a metric space when endowed with the Carnot-Caratheodory distance. 
When the manifold is analytic and the orthonormal frame can be assigned through m analytic 
vector fields, we say that the sub-Riemannian manifold is analytic. This is the case of the problems 
studied in this article. 

Remark 11 In the problem above the final time T can be free or fixed since the cost is invariant 
by time reparameterization. As a consequence the spaces L^ and W^'^ in (9) can be replaced with 
L°^ and Lip (like in (11)), since we can always reparameterize trajectories in such a way that 
ui{t)'^ + U2{t)'^ = 1 for every r € [0,T]. //■ui(r)^ +112(1")^ = 1 for a.e. r G [0, T] we say that the 
curve is parameterized by sR-arclength. See [6, Section 2.1.1] for more details. 

We now state first-order necessary conditions for our problem. 

Proposition 12 (Pontryagin Maximum Principle for 3D contact problems) In the 3D con- 
tact case, a curve parameterized by sR-arclength is a geodesic if and only if it is the projection of 
a solution of the Hamiltonian system corresponding to the Hamiltonian 

H{q,p) = l{{p,X,{q)f + {p,X2{q)f), qeM, pe T*M, (13) 

lying on the level set H = 1/2. 

This simple form of the Pontryagin Maximum Principle follows from the absence of nontrivial abnor- 
mal extremals in 3D contact geometry, as a consequence of the condition dim(Span{Xi, X2, [Xi, X2]}(g)) 
3 for every q G M, see [2]. For a general form of the Pontryagin Maximum Principle, see [4]. 

Remark 13 As a consequence of Proposition 12, for 3D contact problems, geodesies and the cor- 
responding controls are always smooth and even analytic if M, Xi , X2 are analytic, as it is the case 
for the problems studied in this article. Analyticity of geodesies in sub-Riemannian geometry holds 
for general analytic sub-Riemannian manifolds having no abnormal extremals. For more details 
about abnormal extremals, see e.g. [3, 19]. 

In general, geodesies are not optimal for all times. Instead, minimizers are geodesies by defini- 
tion. 

For 3D contact problem, we have that local minimizers are geodesies. Indeed, first observe 
that the set of local minimizers is same if we consider the space W^'^ or W^'°° , see Remark 11. 
Observe now that a local minimizer is a solution of the PMP (see [4]), and due to Proposition 12 
a curve is a solution of the PMP if and only if it is a geodesic. See more details in [3]. 



10 



A 3D contact manifold is said to be "complete" if all geodesies are defined for all times. This is 
the case for the problem (Pmec) defined in Section 4 below. In the following, for simplicity 
of notation, we always deal with complete 3D contact manifolds. 

In the following we denote by {q{t),p{t)) = e^^{qo-,Po) the unique solution at time t of the 
Hamiltonian system 

q = dpH, p = -dqH, 

with initial condition {q{0),p{0)) = {qo,Po)- Moreover we denote by tt : T*M — t- M the canonical 
projection {q,p) i— )■ q. 

Definition 14 Let {M, span{Xi,X2}) be a 3D contact manifold and qo € M. Let Aq^ := {po £ 
T* M\ H{qQ,pQ) = 1/2}. We define the exponential map starting from qo as 

Exp^^^ : Aq, xR+^M, Exp^^^ {po, t) = 7r(e*^(go, Po))- (14) 

Next, we recall the definition of cut and conjugate time. 

Definition 15 Let qo £ M and q{.) be a geodesic parameterized by sR-arclength starting from go- 
The cut time for q{.) is Tcutili-)) = sup{t > 0, | Q'(.)|[o,t] is optimal]. The cut locus from qo is the 
set 

Cut(go) = {QiTcutiQi-))) I Qi-) geodesic parameterized by sR-arclength starting from go}- 

Definition 16 Let qo £ M and q(.) be a geodesic parameterized by sR-arclength starting from q^ 
with initial covector pQ. The first conjugate time of q{.) is 

Tconjili)) = min{t > | {po,t) is a critical point of Exp }. 

The conjugate locus from qo is the set Con(go) = {^(^conj (9(-))) I q{-) sR-arclength geodesic from qo}. 

A geodesic loses its local optimality at its first conjugate locus. However a geodesic can lose 
optimality for "global" reasons. Hence we introduce the following: 

Definition 17 Let qo £ M and q{.) be a geodesic parameterized by sR-arclength starting from qo. 
We say that t^ax > is a Maxwell time for q{.) if there exists another geodesic q{.), parameterized 
by sR-arclength starting from qo such that q{tmax) = <i{tmax) 

It is well known that, for a geodesic g(.), the cut time Tcut{(l[-)) is either equal to the first 
conjugate time or to the first Maxwell time, see for instance [2]. Moreover, we have (see again [2]): 

Theorem 18 Let ^ be a geodesic starting from qo and let Tcut <ind Tconj be its cut and conjugate 
times (possibly +oo^. Then 

• ^ cut _; J conj j 

• 'J is globally optimal from t = to Tcut o,nd it is not globally optimal from t = to Tcut + e> 
for every e > 0; 

• ^ is locally optimal from t = to Tconj cmd it is not locally optimal from t = to Tconj + £, 
for every e > 0. 

Remark 19 In 3D contact geometry (and more in general in sub-Riemannian geometry) the ex- 
ponential map is never a local diffeomorphism in a neighborhood of a point. As a consequence, 
spheres are never smooth and both the cut and the conjugate locus from qo are adjacent to the point 
qo itself, i.e. qo is contained in their closure (see [1]). 
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4 Definition and study of (Pmec) 

In this section we introduce the auxihary mechanical problem (Pmec)- The study of solutions of 
such problem is the main tool that we use to prove Theorem 2. 
We first define the mechanical problem (Pmec)- 

(Pmec): Fix ,f > and (x^, yi„, 6'j„), (x/j„,y/j„, 6'/j„) G M^ ^ S^. In the space of L^ controls 
u{-),v{-) : [0,T] —7- R, find the solutions of: 




{x{0),y{0),e{0)) = {Xin,yin,Oin), {x{T),y{T),e{T)) = {Xfin,yfin,0fi 

t{T) \ I cos{9{r)) \ I 

= u{t) sin(0(r)) +v{t) 

f-T 

I yS.'^u^r)'^ + v[t)'^ dr — t- min (here T is free) (15) 

Jo 

This problem (which cannot be interpreted as a problem of reconstruction of planar curves, 
as explained in [8]) has been completely solved in a series of papers by one of the authors (see 
[18, 24, 25]). 

Remark 20 Observe that (Pmec) (as well as (Pprojective) o,nd (Pcurve)/) depend on a parameter 
^ > 0. It is easy to reduce our study to the case ^ = 1. Indeed, consider the problem (Pmec) 
with a fixed ^ > 0, that we call (Pmec)(C)- Given a curve q{.) with cost C^, apply the dilation 
{x,y) — ;■ {tX, ly) to find a curve q{.). This curve has boundary conditions that are dilations of the 
previous boundary conditions, and it satisfies the dynamics for (Pmec)- If one considers now its 
cost Ci for the problem (Pmec)(1); one finds that Ci = C^. Hence, the problem of minimization 
for all (Pmec) is equivalent to the case (Pmec)(1)- The same holds for (Pprojective); (Pcurve); 
with an identical proof. For this reason, we will fix ^ = 1 from now on. 

Remark that (Pmec) is a 3D contact problem. Then, one can use the techniques given in 
Section 3 to compute the minimizers. This is the goal of the next section. 

4.1 Computation of geodesies for (Pmec) 

In this section, we compute the geodesies for (Pmec) with ^ = 1, and prove some properties that 
will be useful in the following. First observe that for (Pmec) there is existence of minimizers for 
every pair {qin,qfin) of initial and final conditions, and minimizers are geodesies. See [18, 24, 25]. 
Moreover, geodesies are analytic, see Remark 13. 

Since (Pmec) is 3D contact, we can apply Proposition 12 to compute geodesies. We recall that 
we have 

M = R^xS\ q={x,y,e), p = {pi,P2,P3), Xi = { sm{6) |, X 
Hence, by Proposition 12, we have 





-f^ = ;^ ((Pi cos 9 + P2 sin 6*)^ + pi) 
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and the Hamiltonian equations are: 



^ ~ dpi 


= h{q,p) cos 9, 


Pi = 


un 

dx 


= 0, 


^ ~ dp2 


= h{q,p) sine, 


P2 = 


dH 

dy 


= 0, 


e = 


dH „ 


h = 


dH 

de 


= - 



h{q,p){—pi sin 9 + P2 cos 9), 

where h{q,p) = picos9 + p2sin9. Notice that this Hamiltonian system is integrable in the 
sense of Liouvihe, since we have enough constant of the motions in involution. Moreover, it can 
be solved easily in terms of elliptic functions. Setting pi = P^ cos Pa, P2 = Pr sin Pa one has 
h{p, q) = Pr cos{9 — Pa) and 9{t) is solution of the pendulum like equation 9 = \P^ sin{2{9 — Pa))- 
Due to invariance by rototranslations, the initial condition on the q variable can be fixed to be 
{xin,yin,9in) = (0,0,0), without loss of generality. The initial condition on the p variable is such 
that H{0) = 1/2. Hence p{0) must belong to the cilinder 

C = {{pi,P2,P3)\pI+pI = 1}. (16) 

In the following we use the notation of [18, 24, 25]. Introduce coordinates {i^,c) on C as follows: 
sin(z//2) = pi cos^ +P2 sin^, cos(z^/2) = — ^3, c = 2{p2 cos 9 — pi sin 9) , (17) 

with v = 29 + IT G 25"^. Here 25"^ = M/(47rZ) is the double covering of the standard circle 

S^ =M/(27rZ). 

In these coordinates, the Hamiltonian system reads as follows: 

z> = c, c=-sinv, {u,c)e{2Sl)xM.c, (18) 

i; = sin|cos0, y = sin|sin6', = — cos|. (19) 

Note that the curvature of the curve (x(.),y(.)) is equal to 

^=(J^-^ = — tW2)- (20) 

We now define cusps for geodesies of (Pmec)- Recall that both the geodesies and the corresponding 
controls are analytic. 

Definition 21 Let q{.) = {x{-),y{-),9{-)) be a geodesic o/(Pmec); parameterized by sR-arclength. 
We say that Tcusp is a cusp time for q(.) (and q(Tcusp) a cusp point) if u{.) changes its sign at 
Tcusp- We say that the restriction of q{-) to an interval [0,T] has no internal cusps if no t g]0,T[ 
is a cusp time. 

Given a curve q{.) with a cusp point at Tcusp, we have that its projection on the plane x{-),y{-) 
has a planar cusp at Tcusp as well, see Figure 4. More precisely, we have the following lemma. 

Lemma 22 A geodesic j (without angular cusps) has a cusp atTcusp if and only ifliuit^Tcua 1-^(01 - 

00. 
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Proof. First observe that 7 has an internal cusp at Tcusp if) for t — )• Tcusp, it holds u{t) — )■ and 
v{t) 7^ 0, i.e. using (19) one has u{t) = sin(i//2) — ;■ and v{t) = — cos(z^/2) -/^ 0. This is equivalent 
to K{t) = - cot{iy/2) -^ 00, by using (20). D 

Also observe that one can recover inflection points of the planar curve x(-),y(-) from the ex- 
pression of q{.). Indeed, at an inflection point of the planar curve, we have that the corresponding 
q{.) satisfies K = and z/ = vr + 27rn, with n E Z. 

4.2 Qualitative form of the geodesies 

Equation (18) is the pendulum equation 

P = -sinz/, ue2S^ =R/{4TrZ), (21) 

whose phase portrait is shown in Figure 5. 

C 




Figure 5: Phase portraits of the pendulum equation, with the 5 types of trajectories. 

There exist 5 types of geodesies corresponding the different pendulum trajectories. 

1. Type S: stable equilibrium of the pendulum: v = Q. For the corresponding planar trajectory, 
in this case we have {x{t),y{t)) = (0,0). These are the only geodesies with angular cusps. 

2. Type U: unstable equilibria of the pendulum: z^ = vr or z/ = — vr. For the corresponding 
planar trajectory, in this case we have {x{t),y{t)) = (t,0) or {x{t),y{t)) = (— t, 0), i.e. we get 
a straight line. 

3. Type R: rotating pendulum. For the corresponding planar trajectory, in this case we have 
that {x{t),y{t)) has infinite number of cusps and no inflection points (Fig. 6). Note that in 
this case ^ is a monotone function. 

4. Type O: oscillating pendulum. For the corresponding planar trajectory, in this case we have 
that {x{t),y{t)) has infinite number of cusps and infinite number of inflection points (Fig. 7). 
Observe that between two cusps we have an inflection point, and between two inflection points 
we have a cusp. 

5. Type Sep: separating trajectory of the pendulum. For the corresponding planar trajectory, 
in this case we have that {x{t),y{t)) has one cusps and no inflection points (Fig. 8). 
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Figure 6: Trajectory of type R. Figure 7: Trajectory of type O. Figure 8: Trajectory of type Sep. 

The explicit expression of geodesies in terms of elliptic functions are recalled in Appendix A. 

Recall that, for trajectories of type R, O and Sep, the cusp occurs whenever zv(t) = 27rn, with 
n G Z, since in this case one has from Lemma 22 that K{t) — t- cx) for t ^ T^ 



cusp- 



4.3 Optimality of geodesies 

Let q{.) = {x{.),y{.),6{.)) be a geodesic parameterized by sub-Riemannian arclength t G [0, T]. 
Consider the following two mappings of geodesies^: 



S,T: q{ 



i-> 



5(-),9t(-)> 



with q{.) : [0,T] 



xS' 



where 



and 






Xj{t) 

yiit) 



eiT)-9{T-t), 

- COS 6l(T)(x(T) -x{T 

-sme{T){x{T)-x{T- 



- 1)) - sme{T){y{T) 
t)) + cos e{T){y{T) 



y{T 
y{T 



t)), 



0{T-t)-e{T), 

COS 0(T)(x(T -t)- x{T)) + sme{T){y{T - t) - y{T)), 

-sme{T){x{T -t)- x{T)) + cos 0(T)(y(T - t) - y{T)). 



Modulo rotations of the plane (x, y), the mapping S acts as reflection of the curve {x{.)^ y{-)) iii the 
middle perpendicular to the segment that connects the points (x(0),|/(0)) and {x{T),y{T)); the 
mapping T acts as reflection in the midpoint of this segment. See Figures 9 and 10. 

A point q{t) of a trajectory q{.) is called a Maxwell point corresponding to the reflection S if 
q{t) = qsit) and q{-) ^ qsi')- The same definition can be given for T. Examples of Maxwell points 
for the reflections § and T are shown at Figures 11 and 12. 

The following theorem proved in [18, 24, 25] describes optimality of geodesies. 



^Such mappings are denoted by £^,e^ in [18, 24, 25], respectively. 
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Figure 9: Action of S on 1 1— )■ (x(t), yit)). Figure 10: Action of T on 1 1— t- {x{t), y{t)). 





Figure 11: Maxwell point for reflection S. Figure 12: Maxwell point for reflection T. 
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Theorem 23 A geodesic q{.) on the interval [0,T], is optimal if and only if each point q{t), t G 
(0,r), is neither a Maxwell points corresponding to S or T, nor the limit of a sequence of Maxwell 
points. 

Notice that if a point q{t) is a limit of Maxwell points then it is a Maxwell point or a conjugate 
point. 

Denote by Tp(,nd the period of motion of the pendulum (21). It was proved in [25] that the cut 
time satisfies the following: 

• Tcut = \Tpend for geodesies of type R, 

• Tcut G ilTpend,Tpend) for gcodesics of type O, 

• Tcut = +00 = Tpcnd for geodesies of types S, U and Sep. 

Corollary 24 Let q(.) be a geodesic. Let Tcusp, o-nd Tcut be the first cusp time and the cut time 

(possibly +Oo). Then Tcusp < Tcut- 

Proof. For geodesies of types R and O it follows from the phase portrait of pendulum (21) that 
there exists t G (0, \Tpcnd) such that v{t) = 27rn. This implies that K{t) — )■ +oo, and, by Lemma 
22, we have a cusp point for such t. Then Tcusp ^ 2'^pend ^ Tcut- 

For geodesies of types S, U and Sep, the inequality Tcusp < Tcut is obvious since Tcut = +00. □ 

Corollary 25 Let g(.) defined on [0,T] be a minimizer having an internal cusp. Then any other 
minimizer between q{0) and q{T) has an internal cusp. 

Proof. It was proved in [25] that for any points qQ,qi ^ ^ x S^, there exist either one or two 
minimizers connecting go to qi. Moreover, if there are two such minimizers q{.) and q{-), then q{-) 
is obtained from q{.) by a reflection S or T. So if g(.) has an internal cusp, then q{-) has an internal 
cusp as well. D 

5 Equivalence of problems 

In this section, we state precisely the connections between minimizers of problems (Pcurve), (Pprojective) 
and (Pmec) defined above. The problems are recalled in Table 1 for the reader's convenience. We 
also prove that there exists minimizers of (Pcurve) that are absolutely continuous but not Lispchitz. 
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Notation 





X, = I sin 61 , X- 




here x := {x,y) G M^ and 6 ^ S^ or P^ as specified below. We denote with s the plane-arclength 
parameter and with t the sR-arclength parameter. In all problems written below we have the following 

• initial and final conditions ixin,yin,Oin), {x fin, U fin, fin) are given; 

• the final time T (or length i) is free. 



Problem (Pcurve): 

qeR^ xS^ q = Xi+ vX2, / Ve^ + v^ ds = ^^2 j_ K^s^ds -^ 

Jo Jo 

Problem (Pmec): 



mm 



qem^ xS^ q = uXi+ vX2, I -J i^u^ + v'^ dt 

Jo 



mm 



Problem (Pprojective): 

qeR'^xP^ q = uXi+vX2, / y/^^u'^ + v'^ dt = / \\±\\ ^/ ^'^ + K {t^ dt -^ min 

Jo Jo 



Table 1. The different problems we study in the paper. 

First notice that the problem (Pmec) admits a solution, as shown in [18, 24, 25]. The same 
arguments apply to (Pprojective), for which existence of a solution is verified as well, see [6]. 

Also recall that the definitions of (Pprojective) and (Pmec) are very similar, with the only 
difference that 6 £ P^ or 9 £ S^, respectively. This is based on the fact that M? x S'^ is a double 
covering of M? x P^ . Moreover, both the dynamics and the infinitesimal cost in (Pmec) are 
compatible with the projection M^ x S*^ — )• M^ x P^. Thus, the geodesies for (Pprojective) are the 
projection of the geodesies for (Pmec)- Then, locally the two problems are equivalent. If we look 
for the minimizer for (Pprojective) from {xin,yin,Oin) to (x fin, y fin, d fin), then it is the shortest 
minimizer between the minimizing geodesies for (Pmec): 

minimizing geodesic qi{.) : connecting {xin,yin,Oin) to (x fin, y fin, fin); 

minimizing geodesic g2(-) : connecting (xj„,|/i„,6'j„ + vr) to {xfin,yfin,Ofin); 

minimizing geodesic g3(.) : connecting {xin,yin,Oin) to {xfin,yfin,Ofin + T^); 

minimizing geodesic qi{.) : connecting (xj„,|/j„,6'j„ + vr) to ixfin,yfin,Ofin + t^); 

In reality, these four minizing geodesies are coupled two by two: indeed, qi and 94 are geometri- 
cally the same curve, as well as (72 and q^. This is a direct consequence of the fact that one can 
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reparametrize a geodesic backward in time, and as a consequence boudary conditions are trans- 
formed from 6 to 6+it. More precisely, there exists the fohowing symmetry of geodesies for (Pmec): 
{x{t),y{t),e{t)) ^ {x{t),y{t),0{t) +7r), by replacing {u{t),v{t)) with {-u{t),v{t)). See Figure 13. 




Figure 13: Minimizing geodesic for (Pprojective) from minimizing geodesies for (Pmec)- 

It is also easy to prove that a minimizer of (Pmec) without cusps is also a minimizer of (Pcurve)- 
Indeed, take a minimizer q{.) of (Pmec) without cusps, thus with ||x(r)|| > for r G [0,T]. Then, 
reparametrize the time to have a spatial arclength parametrization, i.e. u = ||x|| = 1 (this is 
possible exactly because it has no cusps). This new parametrization of q{.) satisfies the dynamics 
for (Pcurve) and the boundary conditions. Assume now by contradiction that there exists a curve 
q{.) satisfying the dynamics for (Pcurve) and the boundary conditions with a cost that is smaller 
that the cost for q{.). Then q{.) also satisfies the dynamics for (Pmec) and boundary conditions, 
with a smaller cost, hence q{.) is not a minimizer. Contradiction. 

5.1 Connection between curves of (Pcurve) and (Pmec) 

In this section, we study in more detail the connection between curves of (Pcurve) and (Pmec)- 
First of all, observe that (Pcurve) and (Pmec) are defined on the same manifold SE{2). Moreover, 
each curve r(.) = {x{.),y{.),9{.)) satisfying the dynamics for (Pcurve) with a certain control v{.), 
also satisfies the dynammics for (Pmec) with controls u{.) = 1 and v{.). For simplicity of notation, 
we give the following definition. 

Definition 26 LetT{.) = {x{.),y{.),9{.)) be a curve in SE{2) satisfying the dynamics for (Pcurve) 
with a certain control v{.). We define the corresponding curve q{.) for (Pmec) o-s the same 
parametrized curve (x{.),y{.),9{.)), and the corresponding pair as the pair trajectory-control 
((?(.), (n(.),w(.))) withui.) = l. 

We define the corresponding reparametrized pair {qi{.),{ui{.),vi{.))) for (Pmec) the 
time-reparametrization of the corresponding pair {q{-), {u{.),v{.))) by sR-arclength, and the corre- 
sponding reparametrized curve as the curve qi{.)- 

Recall that the time-reparametrization by sR-arclength of an admissible curve for (Pmec) is 
always possible. A detailed explanation for time-reparametrization of a curve with controls in L^ 
to have controls in L°° is given in [6, Section 2.1.1]. 

We now focus on solutions of the Pontryagin Maximum Principle (PMP). For (Pcurve), one 
cannot apply the standard PMP since one cannot guarantee a priori that optimal controls are in 
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L°° . For this reason, we apply a generalized version of the PMP which holds for L^ controls (see 
[28, Thru 8.2.1]). We have the following result. 

Theorem 27 Let T{.) he a solution of the generalized PMP for (Pcurve)- Then the corresponding 
reparametrized curve is a solution of the standard PMP for (Pmec)- 

The proof of this Theorem is given in Appendix B. Here we recall the main steps of the proof: 

STEP 1: we prove that if (r(.), «(.)) is a solution of the generalized Pontryagin Maximum Principle 
for (Pcurve), then, the corresponding pair (g(.), {u{.),v{.))) for (Pmec) is a solution of the 
generalized Pontryagin Maximum Principle. 

STEP 2: we prove that the corresponding reparametrized pair is a solution of the standard PMP. 

We are now ready to discuss the connection between geodesies for (Pcurve) and (Pmec)- 

Proposition 28 Let T{.) be a geodesic for (Pcurve)- Then the corresponding reparametrized curve 
is a geodesic for (Pmec)- 

Proof. Let r(.) be a geodesic for (Pcurve)- By definition, for every sufficiently small interval its 
restriction is a global minimizer. Then it is a solution of the generalized PMP for (Pcurve)- Hence, 
applying the previous Theorem 27, we have that the corresponding reparametrized trajectory is a 
solution of the standard PMP for (Pmec)- This implies that it is a geodesic for (Pmec), due to 
Proposition 12. D 

5.2 (Pcurve) admits minimizers which are absolutely continuous but not Lips- 
chitz 

We now show that the problem (Pcurve) exhibits an interesting phenomenon: there exist absolutely 
continuous minimizers that are not Lipschitz. Other examples are given in [27]. 

Consider a geodesic of (Pmec) defined on [0,T] having no internal cusp and corresponding 
to controls u{-) and f (-). From Corollary 24 it follows that it is optimal. Assume now that this 
geodesic has a cusp at T. Then, by Lemma 22, we have that for t — t- T it holds u{t) — )■ and 
K{t) — )■ cx). Notice that yJu{TY + v{tY is integrable on [0,T], since its integral is exactly the 
Carnot-Caratheodory distance (11), that is finite, see e.g. [18]. Since the cost of (Pmec) and 
(Pcurve) coincide, we have that J„ ^Jl + lC^s)ds is finite. In particular, K{.) is a L^ function that 
is not L°^. Reparametrize time to have an admissible curve r(.) for (Pcurve)) with control v{.). 
Since v{s) coincides with K{s), then v{.) is a L^ function that is not L°^. This means moreover 
that the trajectory r(.) for (Pcurve) has unbounded control and it is not Lipschitz. 

This phenomenon is extremely interesting in optimal control. Indeed, direct application of 
standard techniques for the computation of local minimizers, such as the Pontryagin Maximum 
Principle, would provide local minimizers in the "too small" set of controls L°^([0,T],M). In 
other words, the absolutely continuous minimizers that are not Lipschitz are not detected by the 
Pontryagin Maximum Principle. For this reason, we were obliged to use the generalized PMP for 
(Pcurve) in Theorem 27. 

Instead, the auxiliary problem (Pmec) does not present this phenomenon, since by re-parametrization 
one can always reduce to the set L°^([0, T],M). 
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6 Existence of minimizing curves 

In this section we prove the main results of this paper, proving Theorem 2. We characterize the set 
of boundary conditions for which a solution of (Pcurve) exists. We show that the set of boundary 
conditions for which a solution exists coincides with the set of boundary conditions for which 
a local minimizer exists. Moreover, such set coincides with the set of boundary conditions for 
which a geodesies joining them exists. We also give some properties of such set. 

After this theoretical result, we show explicitly the set of initial and final points for which a 
solution exists, computed numerically. For more details on this subject, see [10]. 

From the following result. Theorem 2 follows. 

Theorem 29 (main result) Fix an initial and a final condition qin = {xin,yin,din) o.nd qfin = 
{x fin, y fin, fin) in M? X S^ . Let q{.) be a minimizer for the problem (Pmec) from q^n to qfin- 
The only two possible cases are: 

1. q{.) has neither internal cusps nor angular cusps. Then q{.) is a solution for (Pcurve) from 

Qin to Qfin ■ 

2. q{.) has at least an internal cusp or an angular cusp. Then (Pcurve) from qin to qfin does 
not admit neither a global nor a local minimum nor a geodesic. 

Proof. We use the notation q{.) to denote trajectories for (Pmec), and r(.) for trajectories for 
(Pcurve)- Recall the results of Section 5. Given a r(.) = {x{.),y{.),6{.)) trajectory of (Pcurve), 
this gives naturally a q{.) = {x{.),y[.),6{.)) trajectory of (Pmec)- On the converse, a q{.) = 
{x{.),y{.),9{.)) trajectory of (Pmec) without cusps gives naturally a r(.) trajectory of (Pcurve), 
after reparametrization. 

Fix an initial and a final condition qin = {xin,yin,Oin) and qfin = (x fin, y fin, & fin)- Take a 
solution q{.) of (Pmec)- If ?(-) has no cusps, then one can reparametrize time to have a curve r(.) 
solution of (Pcurve)- If ?(•) has cusps at boundaries, then the same re-parametrization (that can 
be applied, as explained in Section 5.1) gives the corresponding r(.), that is a solution of (Pcurve)- 
The first part is now proved. 

We prove the second part by contradiction. If q{.) has an internal cusp, then any other solution 
of (Pmec) from qin to qfin has an internal cusp, as proved in Corollary 25. By contradiction, 
assume that there exists f (.), either a solution (i.e. a global minimizer) of (Pcurve) from q^n to 
qfin, or a local minimizer, or a geodesic. In the three cases, the corresponding reparametrized curve 
on SE{2) of (Pmec), that we denote by gi(.), has no cusps. 

We first study the case of geodesies. Let r(.) be a geodesic of (Pcurve)- Then qi{.) is a geodesic 
of (Pmec) between the same boundary conditions of r(.), due to Proposition 28. Then, two cases 
are possible: 

• Let ^i(.) be a solution, i.e. a global minimizer, for (Pmec)- Then both q{.) and q{.) are 
minimizers, one with cusps and the other without cusps. This yelds a contradiction with 
Corollary 25. 

• Let qi{.) be a geodesic for (Pmec) that is not a global minimizer. We denote with [0,T] 
the time-interval of definition of gi(.). Then there exists a cut time tcut < T for qi{.). Then 
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there exists a cusp time tcusp ^ ^cut < T for qi{-), see Corollary 24. Then gi(.) has a cusp. 
Contradiction. 

We have a contradiction in both cases. Thus, if q{.) has a internal cusp, there exists no geodesic of 

(Pcurve) from Qin tO Qfin- 

We now study the case of local minimizers. Let r(.) be a local minimizer for (Pcurve)- 
Then, it is a solution of the generalized Pontryagin Maximum Principle [28, Thm 8.2.1]. Applying 
Theorem 27, we have that the corresponding reparametrized curve ^i(.) is a solution of the standard 
Pontryagin Maximum Principle for (Pmec); and then it is a geodesic by Proposition 12. Since 
r(.) has no cusps, then qi{.) has no cusps either, thus it is a global minimizer. Then both q{.) and 
qi(.) are global minimizers, one with cusps and the other without cusps. This yelds a contradiction 
with Corollary 25. 

Since global minimizers are special cases of local minimizers, we have the result for global 
minimizers too. 

If instead q{.) has an angular cusp, then {xin,yin) = {xfin,yfin), see Remark 3. In this case, 
assume that there exists r(.) either a solution of (Pcurve) (i-e. a global minimizer), or a local 
minimizer, or a geodesic. In the three cases, the corresponding reparametrized trajectory of (Pmec) 
qi{.) must be of Type S, since there are no other geodesies steering gj„ to qjin with {ximVin) = 
(x fin, y fin)- By construction, the solution of (Pcurve) is r(.) = ixin,yin,0{-))- Observing the 
dynamics for (Pcurve) in (2), one has that x,y constant implies that the planar length is ^ = 0, 
then we must have 9in = Of in- □ 

Remark 30 Observe that, as a corollary, we have proved that global minimizers, local minimizers 
and geodesies for (Pcurve) coincide. 

Remark 31 The last part of the proof has its practical interest. It shows the non-existence of a 
solution o/ (Pcurve) in the case of {xin,yin) = {xfin,yfin)- This means that, under this condi- 
tion, it is possible to construct a sequence of planar curves 7"'(.), each steering {xin,yin,0in) to 
ixin,yin,dfin) CL^d such that the sequence of the costs of^^{.) converges to the infimum of the cost, 
but that the limit trajectory 7*(.) is a curve reduced to a point, for which the curvature K is not 
well-defined. See Figure 14- 



6.1 Characterization of the existence set 

In this section, we characterize the set of boundary conditions for which a solution of (Pcurve) 
exists, answering the second part of question Q2. We recall that we just proved that the existence 
set does not change if we consider global or local minimizers or geodesies. 

We prove here some topological properties of such set, and give some related numerical results. 

Proposition 32 Let 5 C M^ x S"^ be the set of final conditions qfin = {x fin, y fin, & fin) for which 
a solution of (Pcurve) cxists, starting from e := (0,0,0). We have that S is closed, arc- connected, 
non-compact. 

Proof. Arc-connectedness and non-compactness of 5 are evident. For arc-connectedness, let qa, qt S 
S. This means that there exist two curves qi{.),q2{-) steering e to qa,qb, respectively. Then the 
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Figure 14: Non-existence of a solution of (Pcurve) for {xin^yin) = (x fin, y fin)- 



concatenation of curves (with reversed time for qi{.)) steers Qa to e to %. For non-compactness, 
observe that ah points on the half-hne (t, 0, 0) are in S. 

We now prove that S is closed. Let x" — )■ x* be a sequence of points in S converging to 



X* G M^ X S"""^. For each x", let (p",t"') be a pair covector-time such that g" 
a geodesic connecting e to x" at time t" (i.e. 
t" < r",„„, where T",„„ is the cusp time for q^{-] 



■■= Expg(p",.) is 
x" = Expg(p",t")) and without internal cusp, i.e. 
c„gp, „..^^^ -cMsp "" """^ v.^^^ U....V. ^^^ y ,^.;. Such pair (p", t") exists because x" G 5. The idea 
of the proof is to pass to a limit pair (p* , t* ) that gives a geodesic without internal cusp steering 
e to X*. The main problem is to prove that p" is a bounded sequence in the unbounded space 
Ag, so that a limit p* exists (eventually passing to a subsequence). This will be a consequence of 
optimality of g"(.) up to t". 

Observe that, due to Theorem 29, since (?"(.) has no internal cusps, then it is globally optimal, 
and in particular the sub-Riemannian distance d given in Definition 9 on the whole M? x S^ satisfies 
(i(e, x") = t". Since the distance is continuous, then d(e,x*) = lim„_>oo ^(e, x") = lim„_^ooi"- In 
particular, t" is a converging sequence. If such limit is 0, then d{e,x*) = implies x* = e £ S. 
If instead the limit is finite but non-zero, we prove that the sequence of the corresponding p" is 
bounded. For each p", consider its three coordinates Pi,P2,P3 with respect to the basis described 
in Section 4.1. Since (p")^ + (^3)^ = 1 in Ag, then p" is unbounded if and only if P2 i^- We prove 
that P2 is bounded by contradiction. Take a subsequence pg ~^ °^' using coordinates (i^",c"') 
for Ag recalled in (17), we have that sin(z/"/2) bounded for all n G A^ implies sin(^"') — )• 0, that 
in turn implies c^ unbounded. By passing to a subsequence, we have c" — t- ±00. Without loss 
of generality, we choose +cxd. In particular, for each A^ > 0, one can find M such that, for all 
n > A^ it holds c"'(0) > M. Take now the dynamics z>" = c'^^c^ = — sin(zv"') given in (18): Since 
|c"| < 1, then t'" > M - t, then z^"(t) - zy*(0) > Mt - t'^/2. Recall that i/" = 2Trk on cusp points, 
and denote with T^^ the first cusp point for the trajectory q^{.)- With a simple estimation, for 



M > 2Wtt one has that there exists t G 



0, M - VM^ - in 



such that u'' 



there exists a cusp time T^^sp ^ -^ 
thus M - \/M2 - 47r -> 0. 



VM^ — Att. Observe now that for A^ 



Since t" < T^i^^p because of optimality, we have t 



(t) - i/*(0) > 27r, then 

-)■ 00 one has M — )■ 00, 
0. Contradiction. 
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Thus p" is a bounded sequence. Passing to a subsequence, we have that (p",t") — )■ {p*,t*). By 
continuous dependence of the solution of the Haniiltonian systems with respect to the initial datum 
p and the time t, we have that x* = Expg(p*,t*). 

We are left to prove that x* = Exp^{p* ,t*) £ S, in particular proving that the cusp time 



The convergence of p" also implies the convergence of 



T*usp for Expg(p*,.) satisfies t* < T*^^^. 

u^-jC^, since they are coordinates^ on the cylinder (16). Without loss of generality, we assume 
ly G [0, 2it]. Since a trajectory has a cusp when i^ reaches 2Trk and <?"■(.) has no internal cusps, then 
i^"((0,t")) C (0,27r), hence iy*{[0,t*]) C [0,27r]. If z^*((0, t*)) C (0,27r), then q*{.) has no internal 
cusps, thus X* G S. We now prove that the other possibility does not occur, by contradiction. 
Assume that there exists t* < t* such that q*{.) has a cusp in t* . Linearize the dynamics (18) 
around i* , that gives = c,c = —v and apply the implicit function theorem: this gives a sequence 
t" —7- i* such that u^it"^) = or 27r, hence q^ has a cusp point in t". Since i* < t* , then for n 
sufficiently large we have t" < t", i.e. q"'{-) is a trajectory with an internal cusp. Contradiction. 

n 

Other properties of 5 (which are evident numerically^) are the following: 



1. all points of S satisfy Xfm > 0; 

2. if qfin G S satisfies 9 fin = vr, then it also satisfies x/j„ = 0; similarly, if qf 
Xfin = 0, then it also satisfies 9 fin = vr. The solutions of a problem with qf 
have a cusp in qfin- 



m ^ S satisfies 
= (0,y/m,7r) 



Figure 15: Geodesies reaching x = 0, upper plane. 



Remark 33 The characterization of S is, in some sense, the continuation of the main results of 
the authors in [6]. There, we proved that there exist boundary conditions such that (Pcurve) did 
not admit a minimizer, i.e. that S is not the whole space SE{2). Here we have described in bigger 



More precisely, since v is an angle, we can always choose a sequence u'^ that is converging. 
® Formal proofs are given in [10]. 
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detail the set of boundary conditions S such that (Pcurve) admits a minimizer, together with proving 
that, given boundary conditions, the existence of a minimizer is equivalent to the existence of a local 
minimizer or a geodesic. 
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A Explicit expression of geodesies in terms of elliptic functions 

In this section, we recall the explicit expressions of the geodesies for (Pmec)- They were first 
computed in [18]. 

The geodesies are expressed in sub-Riemannian arc-length t, and they are written in terms of 
Jacobian functions en, sn, dn, E. For more details, see e.g. [29]. Here {i^,c) are the variables for 
the pendulum equation (18) and (ip, k) are the corresponding action-angle coordinates that rectify 
its fiow: (p = 1, k = 0. See detailed explanations in [18, Sec. 4]. 

Since (Pmec) is invariant via rototranslations, we give geodesies starting from (0,0,0) only. 

Recall that we have classified geodesies of (Pmec) via the classification of trajectories of the 
pendulum Eq. (21), see Section 4.2. We have the following 5 cases. 

• The geodesic of type S has the simple expression q{t) = (0, 0, t). The projection on the plane 
gives the line reduced to the point (0, 0). 

• The geodesic of type U has the simple expression q{t) = (t, 0, 0). The projection on the plane 
is the straight half- line (t, 0). 

• Geodesies of type R have the following expression : 

cos 6{t) = en c^cn (99 + t) + sn (psn (99 + t) , 

sin6{t) = sgn (cos(i//2))(sn(/9cn {(f + t) — cinpsn {if + t)), 

sgn (cos(i^/2)) , ,, , , ,, , ^, , ^, 

x{t) = ^ ,^^ ^^ cn(^(dn99-dn(99 + t)) + snv7(t + E(v?) -Yjiip + t))], 
k 

y{t) = (l//c)[snv?(dnv9-dn(v9 + t))-cnv?(t + E((/7)-E((/7 + t))]. 

• Geodesies of type O have the following expression : 

cos d{t) = A; sn [ip/k)sn {ip + t)/k + dn {ip/k)dn {ip + t)/k, 

sin^(t) = A;(sn {(p/k)dn (99 + t)/k — dn {ip/k)sn (99 + t)/k), 

x{t) = sgn (c)fc[dn {ip/k){cn (ip/k) - en (99 + t)/k) + sn {if/k){t/k + E {if/k) - E ((99 + t)/k)], 

y{t) = sgn {c)[k^sn {Lp/k){cn {if/k) - en (99 + t)/k) - dn {ip/k){t/k + E (99/fc) - E (99 + t)/k)]. 

25 



• Geodesies of type Sep have the foUowing expression : 

cos6{t) = l/(cosh(^cosh((/3 + t)) + tanh(/3tanh((/9 + t), 

sin^(t) = sgn (cos(j^/2))(tanh(y9/cosh((^ + ~ tanh(99 + t)/ coship), 

x{t) = sgn (cos(z^/2))sgn (c) [(1/ cosh 99) (1/ cosh c/? — l/cosh(c^ + t)) + tanh99(t + tanhc/? — tanh(99 + i))], 

y{t) = sgn (c)[tanh 99(1/ cosh 99 — l/cosh(99 + t)) — (1/ cosh 99) (t + tanhc^ — tanh((/? + t))]. 

Pictures of geodesies of type R, O, Sep are given in Figures 6, 7 and 8, respectively. 

B Proof of Theorem 27 

In this appendix, we prove Theorem 27. The structure of the proof is given in Section 5.1. We are 
left to prove STEP 1 and STEP 2. 

STEP 1: If (r(.), v{.)) is a solution of the generalized Pontryagin Maximum Principle for (Pcurve); 
then, the corresponding pair {q{-), (ui-), v{-))) is a solution of the generalized Pontryagin Maximum 
Principle for (Pmec)- 

Proof. Without loss of generality, we provide the proof for ^ = 1 . 

Apply the generalized PMP both to problems (Pcurve) and (Pmec)- For (Pmec), the unmax- 
imised Hamiltonian is T-Lm '■= Pi cos{9)u+p2 sm{9)u+psv + XVu'^ + ^^- For (Pcurve), replace u with 
1: we denote such Hamiltonian with Tic- We denote the maximised Hamiltonians with Hm,Hc, 
respectively. Recall that we study free time problems, thus both the maximised Hamiltonians 
satisfy Hm = and He = 0, see [4, Sec. 12.3]. 

We observe that for both problems there are no strictly abnormal extremals (i.e. solutions with 
A = 0). Indeed, for (Pcurve) abnormal extremals are straight lines, that can be realized as normal 
extremals too. The same holds for (Pmec)- Thus we fix from now on A = — 1 without loss of 
generality. 

Let now {T{.),p{.),v{.)) be a trajectory vector-covector-control satisfying the generalized PMP 
for (Pcurve)- We prove that the corresponding trajectory vector-covector-controls iq{-),p{-), {l,v{.))) 
satisfies the generalized PMP for (Pmec)- The main point here is that T-Lm depends on two param- 
eters {u,v), while Tic depends on v only. Thus, to maximise Hamiltonians, one has more degrees 
of freedom for T-Lm than for Tic- We need to prove that such additional degree of freedom u does 
not improve maximisation of the Hamiltonian. 

We first prove that, if 'v(.) maximises^*^ ^c(9(-),p(-),^(-)), then the choice u{.) = l,v{.) = v(.) 
maximises the Hamiltonian 'Hm{<i{-)-,p{-)-, {u{-),v{.))). First observe that both T-Lm and lie are C°^ 
(except for T-Lm in (0, 0)), and concave with respect to variables u, v and v, respectively. Moreover, 
we have no constraints on the controls. Thus, maximisation of the Hamiltonian is equivalent to 
have VuH = 0. 

We are reduced to prove that q^ = q^ = when evaluated in {q{.),p{.), (l,-i;(.))). Observe 

that ^^ = ^ for u = I; thus, since v{.) maximises Tic, then ^ = 0. Hence ^^ = 0. A 

simple computation also shows that q^ evaluated in u{.) = l,v{.) = v{.) is pi cos{9) +p2sm{6) — 

, _^ , whose expression coincides with He when replacing p3 with its expression with respect to 



^i.c, it maximises T-Lc along the trajectory (g(.),p(.)). 
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the optimal control, that is ps = , " _^ . Since He = 0, then —q^ = 0, hence Hm is maximised by 
ui.)^l,vi.) = vi.). 

Thus we have that Hm = -f^c on this trajectory. Then, since He = 0, then it clearly holds 
TiM {<!{■), Pi-), i^i^i-))) = ^^^ it is also clear that {q{.),p{.)) is a solution of the Hamiltonian 
system with Hamiltonian Hm- Then, {q{-),p{-)) is a solution of the generalized PMP for (Pmec)- 
D 

STEP 2: Let {q{.), {u{.),v{.))) with u{.) = 1 be a solution of the generalized Pontryagin Maximum 
Principle for (Pmec)- Then, the curve reparametrized by sR-arclength is a solution of the standard 
Pontryagin Maximum Principle. 

Proof. Recall that for (Pmec) one can always reparametrize curves by sR-arclength. This also 
transforms trajectoires with L^ controls in trajectories with L°^ controls without changing the cost, 
as explained in Remark 11. Choose such reparametrization. 

As a consequence, a solution to the generalized PMP can be reparametrized to have controls 
in L°^ . Since the expression of the equations are the same for the standard and generalized PMP, 
then this reparametrized curve is a solution to the standard PMP. D 
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